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We demonstrate theoretically the use of genetic learning 
algorithms to coherently control the dynamics of a Bose- 
Einstein condensate. We consider specifically the situation 
of a condensate in an optical lattice formed by two counter- 
propagating laser beams. The frequency detuning between 
the lasers acts as a control parameter that can be used to 
precisely manipulate the condensate even in the presence of a 
significant mean-field energy. We illustrate this procedure in 
the coherent acceleration of a condensate and in the prepara- 
tion of a superposition of prescribed relative phase. 

PACS numbers: 03.75.Fi, 32.80.Qk 



I. INTRODUCTION 

Atomic Bose-Einstein condensates (EEC) are now gen- 
erated almost routinely in the laboratory, and are result- 
ing in numerous applications in fundamental and applied 
science. A key element of this research is the ability to 
coherently engineer and control the state of the conden- 
sate and its dynamics. This is for example the case in the 
generation of atomic solitons [Q-^, which requires that 
a precise phase be imprinted on the condensate; atomic 
four-wave mixing, which involves the splitting of a con- 
densate into three momentum groups Q; the mixing of 
optical and matter waves §-§, which again requires the 
splitting of a condensate; mode-locked atom lasers [ p^ , 
where a condensate wave function is spatially modulated 
by a periodic optical potential, etc. It can be expected 
that the preparation of condensate states of increasing 
sophistication will be required in future applications, in- 
cluding atom lithography [|ll].^ and atom holography 
El- 

The coherent manipulation and control of quantum 
states has been the subject of considerable work in many 
other areas of physics. For example, it is now possible to 
precisely carve the electronic wave function of atoms or 
to excite vibrational modes of molecules p^-p^, using 
precisely engineered optical pulses. Similar techniques 
have lead to spectacular advances in nonlinear optics. 

A common tool to many of these developments is the 
use of genetic algorithms. These multi-dimensional opti- 
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mization techniques proceed by parametrizing a control 
function in terms of a finite set of coefficients, or "genes" , 
a particular set of genes being called a "chromosome." 
The genetic algorithm operates on a set of chromosomes, 
the "population." Its success in achieving a design goal 
is quantified by a "fitness function," a measure of how 
close the action of a particular chromosome is to the 
desired state. The algorithm proceeds by replacing an 
ill-fitted fraction of the population by new chromosomes, 
the "offsprings" , that result from the mating of two par- 
ent chromosomes according to some set of rules; see e.g. 
Refs. [ p"7|Jl8| for details. In addition to that controlled 
combination of chromosomes, random mutations on sin- 
gle genes prevent the algorithm from getting trapped in 
local extrema. The process is iterated until one chromo- 
some reaches a prescribed fitness value. 

Because of collisions, the dynamics of atomic EEC is 
intrinsically nonlinear. Hence, it is difficult in general to 
precisely predict the effects of control fields on the con- 
densate properties. Indeed, condensates are emerging as 
excellent test systems to study quantum chaos p9| , p0[ . 
This is a clear indication that straightforward analytical 
tools are generally unlikely to be sufficient in the coherent 
manipulation and control of BEC. It is therefore natural 
to turn instead to the use of genetic algorithms. The 
main goal of this paper is to illustrate how they can be 
applied to the design of specific momentum states of con- 
densates. 

We specifically consider two examples, the accelera- 
tion of a condensate, and the preparation of a BEC in 
a coherent superposition of two momentum states with 
a prescribed phase difference. The external control is 
provided by two counterpropagating laser fields of ad- 
justable frequency that provide a time-dependent optical 
lattice interacting with the condensate via Bragg scat- 
tering. This is a natural choice, since Bragg scattering is 
a well-established tool of atom optics: It has been used 
in many applications such as the determination of the 
coherence properties of condensates |2^,^, the imple- 
mentation of Mach-Zehnder interferometers to image the 
condensate phase 23 1, the splitting of condensates pi] , 
and the creation of initial states appropriate for nonlinear 
mixing processes [|j . Optical lattices have also been used 
to investigate physical effects such as atomic Landau- 
Zener tunneling p5| , p6t , Bloch oscillations p7|-p9t, and 
the acceleration of BECs Also, a Josephson junc- 

tion array was experimentally realized with a BEC in an 
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optical standing wave |10|,|31|,^ . 

The paper is organized as follows: Section |l| discusses 
the model, establishes our notation, and introduces the 
"quasi-modes" of the condensate used in the subsequent 
analysis. Section III briefly reviews important aspects 
of genetic algorithms. The main results are presented 
in section [V, which illustrates the usefulness of genetic 



algorithms in condensate acceleration and in the prepa- 
ration of macroscopically separated momentum states of 
prescribed relative phase. Finally, section ^ is a sum- 
mary and conclusion. Appendix |A| gives further details 
of the genetic algorithm. 



II. MODEL 

We consider a condensate consisting of N atoms of 
mass M at temperature T — 0, and placed in a frequency- 
chirped optical lattice formed by two counterpropagating 
laser beams of wave vector and central frequency luq. 
The electric field is then given by 



Eiz,t) = 



So 



(1) 



where 5{t) is a time-dependent frequency difference be- 
tween the two beams. We define the detuning A from 
the atomic resonance as A = — i^o: with uja being 
the atomic transition frequency. In order to avoid spon- 
taneous emission we use far-off resonance light, so that 
|A| ^ r, with r the spontaneous decay rate of the ex- 
cited atomic level. We assume that the coupling of the 
light field from Eq. (Q) to the atomic system is character- 
ized by a Rabi frequency = d£o/h, where d is the dipole 
matrix element of the transition. If the Rabi frequency 
is small compared to the detuning, <C |A|, we can 
adiabatically eliminate the excited atomic level, which 
evolves on a much faster time scale than the lower level. 
Making use of the rotating wave approximation which ne- 
glects terms varying at twice the optical frequency wq, the 
resulting time-dependent optical potential for the lower 
atomic level is then given by 



V{z, t) = Vb cospfciz - 5{t)t], 



(2) 



where Vb = ?i|17p/2A is the lattice depth (see also Ref. 
psf for details) . In order not to violate the adiabatic ap- 
proximation, the detuning 5{t) must remain small com- 
pared to the detuning A, and vary slow on the fast time 
scale: 



1^1 « |A| 



« A^. 



(3) 



The instantaneous phase velocity of the lattice fringes 
in Eq. (i. 



(4) 



must remain much smaller than the speed of light at all 
times in order to neglect any time-dependence of the 
wavenumber fci,, as discussed in detail in Ref. [Q. 

The control of the state of the BEC is achieved by im- 
posing a time dependence on the detuning 5{t) between 
the counterpropagating laser beams. It is this time de- 
pendence that is to be determined by the genetic algo- 
rithm. 

We assume that the condensate is tightly confined 
in the transverse direction, so that the system can 
be described semiclassically by a one-dimensional time- 
dependent Gross-Pitacvskii equation (GPE). The con- 
densate evolution in the optical lattice of Eq. (H) is then 
given by 



in—ip(z,t) 



' 2M dz- 



+ V{z,t) 
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NUomz,t)\^i:{z,t), 
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where i^i^, t) is the condensate wave function, Uq = 
iTrh^a/M, and a, the s-wave scattering length reduced 
to one dimension, is taken to be positive. In the presence 
of a periodic potential, one needs to be aware of the possi- 
ble fragmentation of the condensate via a Mott insulator 
transition. Ref. investigates this problem, predict- 
ing the onset of fragmentation when the lattice potential 
depths are in excess of 5-10 lattice recoil energies. Here 
we consider only situations where the lattice potential is 
shallow enough that such fragmentation does not occur. 
Superfluid effects can likewise be ignored, since the crit- 
ical velocity for typical Bose-Einstein condensates was 
recently experimentally determined to lie in the mm/s 
regime [p6["p8| and as we show shortly, Bragg diffraction 
only populates momentum side modes of the condensate, 
spaced by momenta 2?ifci, corresponding to velocities far 
exceeding the critial velocity. 

Our goal in this paper is to manipulate the momentum 
of the condensate. It is therefore convenient to introduce 
the momentum space condensate wave function 
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/2n 



dz ijj{z, t)e 
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Substituting these into Eq. (^) we obtain the correspond- 
ing coupled difference-differential GPEs 
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dkidk2(t>ik - ki + k2,t)(j){ki,t)(t>*ik2,t). (7) 
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The spatial extent of the condensates that we have 
in mind is large compared to the lattice period Tr/fc^. 
Their initial momentum distribution is therefore much 
narrower than k^, Ak <^ k^. From Eq. (^, we observe 
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that the optical lattice couples states separated in mo- 
mentum by fc = ±2fcL, and hence leads to a momentum 
distribution consisting in general of a "comb" of narrow 
peaks of width Afc. Ground state collisions lead to a 
broadening of these peaks, but for small enough particle 
numbers N, it can be expected that this broadening re- 
mains small compared to 2kL. This suggests that it is 
useful to expand the momentum space condensate wave 
function on a basis of "quasi-modes" described by the 
orthonormal mode functions 



Unik) = ^ {e[i2n 



l)kL]-e[i2n + l)kL]} (8) 



(9) 



where 

Ut) -- 



dk Un{k)4>{k, t) 



(2n+l)kL 



and the step function 0(x) is defined as 



e(x) 



: X <0 

1 : x>0 



dk(j){k,t) (10) 



(11) 



The associated "quasi-mode" populations Pn{t) are ac- 
cordingly 



<.(2n+l)fet 

Pn{t)^\Ut)?= / dk\^{k,t)\\ 
J{2n~l)kL 



(12) 



In the following, we use a genetic algorithm to find a 
time-dependent detuning S{t) leading to predetermined 
values of the probability amplitudes Cn{t). Before pre- 
senting selected results of this study, though, we review 
for completeness some important aspects of this opti- 
mization technique. 



III. GENETIC ALGORITHM 

Genetic algorithms proceed by parametrizing a control 
function in terms of a finite set of genes, a particular set of 
genes being called a chromosome. The genetic algorithm 
operates on a set of chromosomes, the population, whose 
action on the system to be controlled is quantified in 
terms of a fitness function. In the situation at hand, the 
control is achieved by the time-dependent detuning S{t), 
expressed by a truncated Fourier series 

m 

5i{t) — cos{i/ujfjt) + bin sm[vuj]^t), 1 = 1,... ,7V, 



i/=i 



(13) 



where lor — hkj^/M is the recoil frequency. Each de- 
tuning Si is encoded in a chromosome Ci consisting of 



n = 2m genes ^y , each gene corresponding to one par- 
ticular Fourier coefficient, 



Ci{gii, ■ ■ ■ ,5m) = Ci{aii, . 



1 0>r 



,bii, 



(14) 



The index i labels a specific chromosome, and the size of 
the chromosome population is Af. 

Starting from a randomly initialized population, the 
genetic algorithm uses a set of mating rules, mutations, 
and a problem-specific fitness function /(q) to create 
new generations of chromosomes, as illustrated in Fig.^ 
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Evaluate fitness f(c ;) based 
on <!> i C^f ) , rank all c ; 
by increasing fitness 



Select mating operator, 
select parents Ci and cj 
according to ranking. 

Produce new offspring matrix C; 



FIG. 1. Schematics of the genetic algorithm: In our prob- 
lem, the fitness is evaluated by evolving an initial momentum 
space wavefunction (;A(0) according to the GPE from Eq. (P) 
for a time tf, the dynamics of the optical lattice being deter- 
mined by the detuning Si{t). The final wavefunction is then 
compared to the optimization goal. 

The first step consists in selecting parent chromosomes 
that are to be combined by mating operators. This is 
achieved by ranking the initial Af chromosomes according 
to their fitness and using the so-called "roulette wheel" 
method to preferentially select parent chromosomes 
with a high fitness. 

In the next step mating operators are selected that 
generate a group of "offspring" chromosomes from the 
"parent" population a. As such, these operators are at 
the heart of the genetic algorithm. 
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Several mating operators may be considered: The 
"one-point-crossover" operator cuts the two parent chro- 
mosomes at a randomly chosen position fj, and swaps the 
genes according to 

ci(ffii, • ■ ■,gin) 

C2(ff21, • ■ ■,g2n) 
i 

1- point-crossover (15) 

i 

c'ligii, ■ ■ ■ , gifi, 92,1^+1, ■ ■ ■ , g2n) 

C2(521, • ■ • , g2ti, 9l,ti+l, ■ ■ ■ , 9ln), 

with 1 < n < n. A slightly modified version of this 
operator is the "two-point-crossover" operator that cuts 
the two parent chromosomes at two random positions /ii 
and fi2 and then exchanges the genes between these two 
positions: 

Cl(ffll, • ■ • ,5lri) 
C2(ff21, • ■ • ,52n) 

i 

2- point-crossover (16) 

i 

c'i(5iii • • ■ 1 9l^ll: 92,fii+i, ■ ■ ■ , 92412-1, 9iti2-: ■ ■ ■ :9in) 

^2(^21, ■ ■ • ,52^111 51,^11+1: ■ ■ • i5l,P2-l7 52Ai2J • • • :92n), 

with 1 < fii < fj,2 < n. 

Another type of mating operator, the "average- 
crossover" operator, produces just one offspring from the 
two parent chromosomes by averaging the genes between 
two randomly chosen positions /ii and fi2'- 

ciigii, ■ ■ • ,5ln) 

C2(ff21, • ■ • ,52n) 

i 

average-crossover (17) 
i 

Cl(5ll: ■ • ■ 5^11 + 1: ■ • ■ i 5^i2 - 1 ' S'Ia'2 i ■ • ■ i9ln), 

with 1 < Ml < A*2 < ?i and g'^ = (gi„ + g2K)/2. 

Except for the random location at which the splicing of 
the chromosome occurs, the mating algorithms discussed 
so far are deterministic. In addition, genetic algorithms 
also require the use of "mutation" and "creep." These 
are random operators that produce one offspring from 
one parent by altering any gene of the parent chromo- 
some with a given probability called the mutation rate, 
respectively the creep rate: 

c{gi, ■ ■ • ,5«) 
i 

mutation, creep (18) 

i 

c'{gi,---,g'f,,---,gn)- 



The mutation operator chooses g'^ randomly within some 
bounds, whereas the creep operator shifts the old value 
by a random amount, g'^ = -|- (0.5 — r)pcrocp- Here 
Pcrccp is a parameter that controls the range of the shift 
and < r < 1 a random number. 

The operators discussed in this section are the only 
ones used in our analysis. When it comes to selecting a 
specific mating operator, there are basically two possi- 
bilities: The first one consists of assigning fixed weights 
to the available operators and then choosing randomly 
among them. This is a straightforward approach, but 
it suffers from the problem of not discriminating against 
mating operators that do not perform well for the op- 
timization problem at hand. Thus, a second possibil- 
ity is to dynamically adjust the operator weights over 
the course of the optimization [HTpfl]. This guarantees 
that the best suited operators are applied and allows one 
to test the performance of new mating operators. This 
is done by assigning an adjustable "operator fitness" to 
each mating operator under consideration. As such, the 
mating operators are selected according to their fitness 
the same way as the parent chromosomes are picked. The 
details of the procedure used in our simulations are dis- 
cussed in the next section and also in the appendix. 

IV. RESULTS 
A. Coherent acceleration 

We now apply the genetic algorithm approach to the 
manipulation of the state of a BEC in a chirped opti- 
cal lattice. As a first example, we consider the coherent 
transfer of a condensate population between the adjacent 
quasi-modes uq and ui. Our motivation here is the need 
to find efficient ways to accelerate condensates for atom 
optics applications. This problem was theoretically an- 
alyzed in Refs. and experimentally demonstrated 
in Ref. for the case of a linear dependence of 6{t), 
6{t) — rjt. While this chirping of the lattice detuning 
does lead to a large mean acceleration of the condensate, 
it unfortunately leaves a substantial fraction of the con- 
densate in lower quasi-modes. This translates into a loss 
in coherence that is unacceptable, e.g., in interferometric 
applications. The question, then, is whether an opti- 
mized time-dependent detuning determined by a genetic 
algorithm can eliminate this problem. 

For the sake of illustration, we assume that the con- 
densate is initially in the zero-momentum quasi-mode, 
Po(i = 0) = 1, and seek a time-dependent detuning such 
that Pi{tf) = 1 after some predetermined time tf. In 
that case, the algorithm fitness has the simple form 

/(c,)=Pi(t/), (19) 

with an optimal value of unity. 

Fig. 1^ summarizes the results of the optimization pro- 
cedure. It compares the optimized population transfer 
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of a fairly large condensate consisting of 10^ atoms to 
the case of the same condensate subject to on-resonance 
Bragg scattering. In this example, the genetic algorithm 
involved a one-point-crossover, a two-point-crossover and 
an average cross-over mating operator. In addition, it in- 
cluded two mutation operators with mutation rates of 
0.8 and 0.4, and two creep operators, both at a rate 
of 0.9, but different creep parameters, a "coarse" creep 
with Pciocp = 0.01 and a "finer" creep operator with 
Pciocp = 0.001. More details of the simulations are given 
in the appendix. 




t[ns] 

FIG. 2. Momentum transfer of a Sodium condensate with 
a Gaussian spatial distribution of longitudinal width 50 fim 
and transverse width 5 fim; Xl = 985 nm and Vo = huiR, 
M = 3.82 • lO"^** kg and = 4.9 nm. (a) The solid line 
is the optimized time-dependent detuning. The dotted lines 
label the first and second Bragg resonance (BR) for refer- 
ence, (b) Temporal evolution of the mode population pi for 
iV = 2 ■ 10^ atoms: On resonance (dashed) and optimized 
(solid). The vertical dotted line denotes the time it takes to 
resonantly transfer all population to the quasi-mode Ui in a 
linear two-mode system. 

While it can be expected that resonant Bragg scatter- 
ing at the appropriate frequency transfers perfectly the 
population of a small condensate from mode uq to mode 
Ml, such is not the case for the large condensate we inves- 
tigated (dashed line in Fig.^). In this case, the mean-field 
nonlinearity of the condensate is no longer negligible. It 
dynamically shifts the Bragg resonance so that the 
transfer efficiency drops to barely over 90% and the max- 
imum transfer occurs later in time. It is in such nontriv- 
ial situations that genetic algorithms are expected to be 
useful. Indeed, the optimal time-dependent detuning d{t) 



found by the genetic algorithm is highly non-trivial. The 
temporal dependence of the detuning 6(t), which trans- 
fers more than 99% of the population to the quasi-mode 
Ml, reveals that while it is initially advantageous to re- 
main close to the Bragg resonance frequency, as indicated 
by the rather flat portion of the detuning, it eventually 
becomes necessary to drastically couple the condensate 
atoms to higher momentum modes so as to drag the re- 
maining population to the final state mi ]4^ . 

The effect of the mean-field energy is further illustrated 
in Fig. ^ which shows the final momentum distribution 
(f>i{k) within the quasi-mode mi for various numbers of 
atoms in the condensate for optimal transfer. While this 
distribution remains extremely narrow compared to the 
quasi-mode width 2kL, collisions lead to a substantial 
reshaping and broadening within that mode. 




k [units of kj^] 



FIG. 3. Effect of the nonlinearity on the momentum space 
wave function: Momentum space densities within quasimode 
iti after optimizing the transfer, so that pi > 0.99 after time 
tf = 43 us; N = 1-10^ (solid), N = 5-10^ (dashed), iV = I-IO'^ 
atoms (dotted), and A*' = 2 • lO'' (dashed-dotted). 



B. Coherent superposition 

In a second application, we set out to design an equal- 
weight superposition of the two quasi-mode states uq and 

Ml 

0(fc, t/) = (Moe''-^" -I- Mie^'^i) , (20) 

with a prescribed relative phase Atp = cpi — i/jq- In con- 
trast to the previous example, we now want to control 
two properties of the quantum state, the relative phase 
as well as the population in each state. The fitness func- 
tion to be optimized is therefore more complicated. 

We choose to employ a so-called "penalty function" 
P{ci) for the optimization of the quasi-mode popula- 
tions jl^. The goal of P{ci) is to decrease the fitness of 
chromosomes that do not fulfill the desired requirements, 
thereby steering the population towards the target val- 
ues. A prototype penalty function is 




Po{ci),pi{ci) > 0.465 
Po{ct),pi{c,) > 0.47 
Po{c^),pi{c,) > 0.475 
else 



(21) 



The fitness function for tliis optimization problem is then 
given by 



/(cO = 1 - |a(c,) - Aip\ - P{ci), 



(22) 



where, as we recall, Ci corresponds to a specific realiza- 
tion of the time-dependent detuning 5(t) and a(ci) is the 
relative phase corresponding to this realization. This fit- 
ness function reaches its maximum, unity in this case, 
when the populations are within the specified range and 
the phase difference is exactly as prescribed. The results 
for the optimization for the two cases At/? = — 7r/2 and 
Aiy9 = — 7r/4 are shown in Fig. ^. For the genetic algo- 
rithm we used the same operators and their parameters 
as in Sec. |IV A . More details of the simulations are given 
in the appendix. 

Fig. 1^ shows the momentum distributions (/>o(fc) and 
4>i{k) of the condensate within each of the two quasi- 
modes uq and ui, as well as the corresponding phases 
for the case A(p — —tt/2. Clearly, the genetic algorithm 
converges to the stated goal, and produces a condensate 
in the desired coherent superposition. 

In particular, we observe that the relative phase of the 
two components is approximately constant in the region 
where the condensate wave function is different from zero. 
The optimization goal of Aip = —tt/2 is achieved with an 
accuracy of over 99% at the center of the mode. The 
curvature of the phase at the wings of the wave function 
is due to nonlinear phase shifts accumulated during the 
Hamiltonian evolution. It could be reduced by decreasing 
the number of atoms in the condensate. 




t[ns] 



FIG. 4. Excitation of an equal-weight coherent superpo- 
sition of momentum states with relative phases, Aip — — 7r/4 
(solid) and Aif = —n/2 (dashed). The parameters are the 
same as in Fig. ^, except that the condensate has a Gaus- 
sian spatial distribution of longitudinal width 100 fim and 
N = 5 ■ 10*. (a) Optimized time-dependent detuning used to 
create the superposition, (b) Temporal evolution of the mode 
population of the two involved modes, (c) Temporal evolu- 
tion of relative phase of the two quasi-modes at the center of 
their momentum distribution, see Fig. 5. 
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k [units of kj^] 
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-0.0066 -0.0033 0.0033 0.0066 
k [units of kj^] 

FIG. 5. Density profiles and phases of the two momentum 
modes after the optimization at t = tf. The lower abscissa 
corresponds to the momentum wavefunction 00 with phase 
ipo centered around fc = 0, the upper to (j}i with phase ipi 
centered around k — 2fci. Note that pi{tf),p2{tf) > 0.475. 



V. SUMMARY AND CONCLUSION 

In summary, we have demonstrated the usefulness of 
genetic algorithms in the control and manipulation of 
the quantum states of Bose-Einstein condensates. This 
was illustrated in two examples, the coherent population 
transfer between momentum states in a large condensate 
where the mean-field effects are important, and the cre- 
ation of coherent superpositions of states of prescribed 
population and relative phase. We found that time- 
dependent Bragg scattering combined with the powerful 
optimization capabilities of genetic algorithms, provides 
a novel tool for quantum state design and coherent con- 
trol in linear and nonlinear atom optics. The extension 
of these ideas to integrated atom optics appears partic- 
ularly promising. Future theoretical work will general- 
ize these concepts to systems with more degrees of free- 
dom, such as e.g. multi-component condensates, and to 
additional control mechanisms such as time-dependent 
magnetic fields. The extension of this work to quantum- 
degenerate Bose-Fermi mixtures also appears promising. 
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APPENDIX A: DETAILS OF THE GENETIC 
ALGORITHM 

In this Appendix we give specific details of the im- 
plementation of the genetic algorithm in the problem of 
BEC state engineering in an optical lattice. 

Over the course of the optimization we monitor the 
maximum fitness /max of a population. 



/max = max{/(ci) : i = 1, . . . ,7V}. 



(Al) 



Obviously, if at least the best chromosome is kept from 
the old generation, the maximum fitness is a monotoni- 
cally increasing function. This feature, called "elitism", 
is used throughout our simulations. Another observable 
of interest is the mean fitness /moan of a population. 



/. 



/ mean 

1 ^ 

mean = ^ ^ f i^i) • 
i—1 



(A2) 



A typical evolution of these two quantities is shown in 
Fig. |^(a) . The maximum fitness increases monotonically 
to reach a value close to the optimum after about 40 
generations. In contrast, the mean fitness rises over the 
course of the first 10 generations and then exhibits fluc- 
tuations due to the stochastic character of the genetic 
algorithm: Parts of the population are replaced by ran- 
domly created chromosomes from one generation to the 
next and randomize the mean fitness value. In our sim- 
ulations we used populations of size J\f = 50 — 100 and 
performed the optimization over 50 — 100 generations. 
We always kept the best chromosomes of a generation 
and replaced 80% — 90% of the population by newly gen- 
erated chromosomes. For the population transfer of Sec. 
EVA we used 16, and for the superposition state engi- 



neering of Sec. [VB 26 genes per chromosome. The gene 
boundaries were chosen as — 0.7u;j^ < a^, bi^ < 0.7wfl,. 

As mentioned in Sec. we use an adaptive oper- 
ator technique, where the operators themselves are dy- 
namically assigned a fitness based on their performance. 
Choosing a particular mating operator via a roulette- 
wheel method Q then assures that good operators are 
employed more often in the mating process. If any op- 
erator produces an offspring that is better than the best 
chromosome of the previous generation, we reward this 
operator by giving it a credit proportional to the increase 
in fitness it caused. Also, we pass half of this given credit 
back to the operator that created the parent chromosome 
involved in producing the better offspring. Thus opera- 
tors that perform well and also their direct ancestors can 
accumulate credit over the run of the simulation. Passing 
credit back to previous operators enables us to reward 
pairs of operators that work well together at a certain 
stage of the optimization process. In our simulations we 
adjust the fitness of all operators every five generations 
based on the credit they accumulated during that period: 
The new operator fitness is then a weighted sum of the 
old fitness (the "basis portion", in our case 85%) and 
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the accumulated credit (the "adaptive portion", in our 
case 15%). The more credit an operator accumulates the 
higher will be its new fitness and the more likely it is 
that it will be chosen in future mating processes. Since 
the total operator fitness is set constant, we introduce a 
lower bound of 0.1 to the operator fitness, thereby pre- 
venting operators that do not perform well over several 
generations from being practically expelled from the pool 
of operators. 

Fig. H(b) shows a typical evolution of the fitness of the 
individual mating operators in t he su perposition state 
engineering problem from Sec. iVB. The creep and 



one-point-crossover operators perform well for the first 
30 generations. For subsequent generations the two- 
point-crossover and the average-crossover operators take 
over and help increasing the maximum fitness, which ap- 
proaches its optimum value (unity in the present exam- 
ple). There is no further improvement in the remain- 
ing 50 generations, the fitness of the various operators 
staying constant. The random mutation operators never 
perform well in the problem at hand. Consequently their 
fitness is quickly reduced to the lower bound. The creep 
operators, which are basically controlled mutations, per- 
form much better and could replace the pure mutation 
operators. 
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FIG. 6. A typical run of the genetic algorithm for 
the superposition state engineering: (a) Maximum fitness 
(solid) and mean fitness (dashed) of the population as a func- 
tion of the generation. We started the optimization with 
pre-optimized chromosomes from previous simulations, which 
explains the high initial maximum fitness of over 60%. (b) 
Operator fitness as a function of the generation, starting with 
equal fitness for all operators. The minimum operator fitness 
is set equal to 0.1 for all operators. 
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